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PREFACE 


This report constitutes the first of two volumes which summarize 
the work accomplished under Contract No. NAS8-21282. It includes sup- 
porting experimental data along with an analysis which is designed to 
predict natural frequencies and forced axisymmetric response in a vehicle 
propellant tank. The results are required for further analysis of dynamic 
stability, which has been accomplished in the second phase of the work 
under the above contract, and published in Final Report, Part II, entitled 
"Dynamic Stability of Cylindrical Propellant Tanks. " The latter document 
includes a listing of a digital computer program developed to facilitate use 
of results from the entire research effort. 

Both Part I and Part II of this Final Report are published on the 
same date. They present significant extensions and refinements of concepts 
originated under previous investigations conducted for NASA-MSFC. Results 
of this preliminary work are summarized in "Dynamic Stability and 
Parametric Resonance in Cylindrical Propellant Tanks, " by Daniel D. Kana, 
Wen-Hwa Chu, and Tom D. Dunham, Final Report, Contract No. NAS8- 
20329, SwRI Project No. 02-1876, January 17, 1968. 
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ABSTRACT 


The influence of a localized rigid mass has been determined for the 
longitudinal response of a model vehicle propellant tank. A theoretical 
analysis is conducted to predict the dynamic response of a simplified model 
consisting of a partially liquid filled cylinder which has a rigid flat bottom, 
an internal ullage pressure, and a rigid top mass. Donnell shell equations, 
along with additional terms to allow for geometrical nonlinearities, are 
used to determine natural frequencies and forced axisymmetric response 
of the model. Numerical results are compared with experimental observa- 
tions for a range of several significant parameters. It is found that strong 
coupling occurs between the motion of the shell, liquid, and top mass 
over a wide portion of this range. 
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PRINCIPAL NOTATION 


a 

a 0 

c 0 


E 

f 

g 


8x 


radius of the shell 
inner radius' of the tank 
speed of sound in the liquid 

speed of stress waves in the shell, (E/p s )l/2 
modulus of elasticity- 
excitation frequency in cps 
standard acceleration of gravity 
noridimensional excitation amplitude , XQto^/g 
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M s 
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Pr 


nondimensional liquid depth, h/a 
nondimensional thickness of shell, h s /a 
depth of liquid 
thickness of shell 

mass moment of inertia of top weight about z axis 
length of the shell 
top mass 

shell mass, 2irah s fp s 
nondimensional top mass 

one -half of the number of circumferential nodes ; cos (m0) 
static stress resultants 
axial wave number, sin mrx/f 

nondimensional pressure loading on cylinder, p r /E 
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r, e, X 
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U,V,W,X,Y 


U, V, w 

x 0 

xo(t) 

x 0 


p 

P s 
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PRINCIPAL NOTATION (Cont'd) 

pressure in the liquid 

ullage pressure 

pressure loading on the shell 

generalized loading due to bottom motion in a rigid tank 

cylindrical coordinates (space fixed) 

time 

u, v, w, x, y, nondime nsionalized by the radius a 
displacements along x, 0, r direction, respectively 
nondimens ional amplitude of axial excitation, Xq = xp/a 
excitation displacement 

displacement amplitude of axial excitation 

axial wavelength parameter , nua/f 

Poisson’s ratio 

mass density of liquid 

mass density of the shell 

nondim en s ional f r equenc y , wa/c s 

k-th eigenvalue, co^a/c s 

circular frequency of excitation 


Superscripts and Subscripts 

( )P related to the forced response 

( A ) the amplitude of ( ) 
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INTRODUCTION 

The interaction of liquids and their elastic containers plays an impor- 
tant role in overall longitudinal response of launch vehicles since liquid pres- 
sures form complex effective masses, as well as avenues of energy feedback 
in a vehicle structure. Asa result, synthesis of systems into relatively 
simple spring-mass representations* provides only a gross approximation of 
vehicle response, although such procedures are very useful from a design 
point of view. 

Various types of responses which occur in a typical single propellant 

2 

tank have been described by Kana and Gormley , while a membrane approach 

* • ^ 
to axisymmetric responses has been given by Kana, et al. Both of these 

articles contain substantial reference lists for additional studies which have 
been reported. All of these previous studies which include a more exact ap- 
proach to the liquid-tank interaction problem have concerned only single pro- 
pellant tank systems, although internal pressurization and elastic support 
effects were included. 

The purpose of the present work is to extend the analysis of a single 
tank to include a rigid top mass, so that a better understanding can be ob- 
tained for the coupled effects which occur in an overall vehicle structure. 

Both symmetric and nonsymmetric responses are considered, although the 
latter are shown to be unaffected by the magnitude of the top mass. Symmetric 
responses are carefully considered, in view of their importance in spring- 
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mass model representations , as well as their formation of the parametric 

4 

loadings for dynamic stability in liquid-tank systems . 

It will be seen that the problem under concern is complicated by the 
fact that the eigenvalues appear in the boundary condition at the top of the 
tank for the case of symmetric responses. In order to handle this cumber- 
some feature, a technique which includes Fourier expansions combined 
with polynomial functions is utilized to produce a matrix-eigenvalue 
formulation. 

THEORETICAL ANALYSIS 


General 

A schematic of the system to be investigated is shown in Figure 1, 
where some parameters are also given for the experimental apparatus. Thus, 
we consider a partially liquid filled thin cylindrical shell, which is internally 
pressurized, has a rigid flat bottom, and supports a rigid top mass. First, 
we consider the case of both symmetric and nonsymmetric free vibration, 
and then axisymmetric response to longitudinal excitation. 

The governing shell equations for the present problem can be obtained 

5 6 

from those developed by Sanders ’ , after retaining only those terms which 

result from geometric nonlinearities. Thus, the effects of static stresses 
will be accounted for. The equations are 

9^u , v 8^v , v 8w , 1 - v / 9^u , 8^v \ _ 1 -_v^ ( n v n „ ^ 

dx 2 + a 8x80 + a 8x + 2a 2 V 80 2 + a 8x90 j " Eh g V shs 9t 2 / ( ] 
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(lc) 


where N xxs , N x q sj Nqq s are axisymmetric static stress resultants formed 
by internal pressurization, hydrostatic liquid pressure, and the static weight 
of the top mass. 

Natural Frequencies of Partially-Filled Tank with Top Mass 

This section deals with vibrational modes whose radial displacements 
are proportional to cos m0, The case of axisymmetric vibration is contained 
as a special case with m = 0. Expanding the axial displacement into a Fourier 

^i> 

cosine series in x plus a quadratic polynomial to satisfy the two end condi- 
tions, and circumferential and radial displacements into Fourier sine series 
in x, there results: 


2 2 

*Note that this procedure is equivalent to expanding 3 U/ 3X into a Fourier 
cosine series. 
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U = [ 1 b 2 X 2 + B, (x - L) + B m0 + f B mn cos (Spjl cos (m«) cos (ut) 

n = 1 ' ' J 


oo 


V= Z c mn sin 1^7^ ) sin (m0) cos (wt) 
n= 1 ' 1 ' 

°° f \ 

w = £ A mn sin ( cos (m0)cos (ut) 

n = 1 \ / 

The boundary conditions are: 

At X * 0: 


0 < x < i 
0 £ 6 S. 2ir 


(2a) 

(2b) 

(2c) 


U - 0, W = 0, V * 0 , 

At X * f /a: 


d 2 W 

ax 2 


= 0 


|~=M**fi£u, W = 0, V = 0 , 


( 3 a, b, c, d) 


(4 a,b, c, d) 


ax‘ 


where 


M(1 - v 2 ) 

M** = — - — - 

2irp s a 3 H s 


M** = 


1 (1 - v2) Ij 


Ps h s a ' 


M** = oo 


for m = 0 

for m = 1 
for m > 2 


(5) 


Equations for determining the unknown coefficients in Equations (2) 
willnowbeformulated. To satisfy the boundary conditions on axial displace- 
ment, one has, from Equation (3a), 


oo 


B 1 a " B m0 " Z B mn - 0 

n= 1 


( 6 ) 
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Boundary condition (4a) requires 

B2 (r )(»5I>nEr- 0 i [( 7 ) B2 (f) 2 + B -o 


+ I B mn (-l) n =0 (7a) 

n = 1 


for m = 0, 1 or 


( 2) B2 (f ) 2 + Bm0 + % B mn(-l) n = G 


for m > 2. 


Upon substitution of Equations (2) into the axial equation of motion 


(la), the x-independent Fourier coefficient (n 1 = 0) yields 


i ? (1 - V] 

1 - m 4, ' * 


Xoa + m 


2 Hull 


— IM0-- -1 


X i a - — 1 B i + m 


2 t 1 - v) 


where 


- £2^(1 ' ^ 2 ) ( B 2 x 20 + B 1 ( x 10 " ”) + B m0 - 0 


- G) (; 


Similarly, for n' > 1, the cos (n'lrx/f ) coefficient of the axial 


Equation (la) yields 


°° / V OO 

Z ^ x n + — 2~ m / Bmn ^ n ' n ” Z 


rnX nCmnSn'n 


00 1 - / 1 \ 
v I X n A mn6n'n + mZ [j B 2 x 2n' + B l X ln' j 

n = 1 

t oo - 

i. X2n' B 2 + Xin'Bi + Z BmnSn'n = 0 
2 n = 1 J 
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The sin (n'trx/i ) coefficient of the circumferential Equation (lb) yields 


oo 


2 m ^ln'®2 + ^On'®l) “ m Z ^ Mfi'n^mn 

n - 1 


oo . oo 

+ Z 2 ^ * m y *“' mn ^ n ' n ^ m Z A mn 6 n 'n 


n = 1 


n = 1 


oo 


- (1 - v ) £2^ Z ^mnSn'n - 0 
n = 1 


(10) 


The sin (ii'ttx/JII ) coefficient of the radial Equation (lc) yields 


, _ go \ °° 

( ^ln'®2 + ^0n'®l " Z ^n®n'n®mn + m Z ^mn^n'n 

V n = 1 ‘ n = 1 


oo 

+ z 

n = 1 


TJL 

, , s . 2 , 2*2 

1 + 7Z (Xn + ’ 


5n'n 


1 - 

Ho 


p 0 2 
2 E Xn 


_l_ 

2 ira^E 


-1 ^n-n 

E J nn 


+ £ ^~m 2 Xo„i ri > A. 


Bn 1 !! f ran 


- (1 - v 2 ) Z A mn( 6 n'n + M mn'n) " 0 

n = 1 


oo 


( 11 ) 


where, for Equations (9) through (11), we define 


- nTra 
x n ~ a 


n T Tra 
x n' _ i 


JL I J 

^In '*ff a X S in( 2 i 2 .x) dX , Son'=T/ sin X ) 

0 0 


dX 


6 J j 

~ J Xcos dX , X2n’ = y - f X 2 cos X^ dX 


Xln' = 
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X 3n'n = T~ / (H " X > sin <X n X > sin <V X > dX 
0 

and, in Equation ( lc), for the added mass pressure of the cosmd mode 
[obtained from Equations (1 8) in the next section] we have used 

oo D oo oo f 

p r = I = I I H s n k A mn M mn'n sin x ) cos (md) 

Thus, in Equation (11), the effect of initial tension due to the ullage 
pressure, the hydrostatic pressure, and top weight are included in addition 
to that of the dynamic liquid pressure. 

Equations (6) through (11) can be cast into the following matrix form: 


[Uj] [Vj] [ Wj ] {o} {y x } j} 


{ A mn} 


[Ril loi [o] {o} {o} {o} 


{^mn} 

[U 2 ] [v 2 l [w 2 | {o} {y 2 } {z^ 


{ B mn} 


[o j fs 2 l f o 1 {o} {p 2 }{q 2 } 


{®mn} 

[U 3 ] lv 3 ] [w 3 ] {o} {y 3 } {z 3 } 


{ C mn} 

■' 


[ 0 ) to 1 l T 3 1 {o} {o} {o} 


{ C mn} 

(0) (V 4 ) (0) x 4 y 4 0 


B mO 

’ - fi k 

(0) (0) (0) 000 

< 

B m0 

(0) (V 5 ) (0) x 5 Y 5 Z 5 


B 1 


(0) (S 5 ) (0) n 5 0 q 5 


B 1 

( 0 ) ( 0 ) ( 0 ) x 6 y 6 z 6 


B 2 | 

N. J 


(0) (0) (0) n 6 p 6 q 6 


b 2 

C * 

3N + 3X 3N + 3 

3N + 3 X 1 

3N + 3 X 3N + 3 

3N+3 X I 


where all series have been truncated to N-terms. To form Equation (12), we 
have written, respectively, Equations (11), (9), (10), (6), (7), and (8) in this 
order. Further, 

hi = (1 - v 2 ) 

while [ ] designates a square N X N submatrix, j | an N X 1 column 
submatrix, ( ) a 1 XN row submatrix, and all other elements are scalars. 
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Note that a [0] indicates a square submatrix whose elements are all zero, 
etc. Elements of all these submatrices are given in the Appendix. 
Submatrices of the 3N + 3X1 column matrices are N X 1 column matrices 
given by 


r A m l ' 


r B m l" 


r r* 

^ml 

A m 2 

, •* 1 

“ > = < 

B m2 

> • =■* 

Cm2 

k A mN j 

1 



, c mN, 


(12a) 


Equation (12) can be expressed as 

[m - ng iri] {a} = o 

Then, since ([r]J is a singular matrix, the eigenvalues and eigenvectors 
can be determined from a standard computer subroutine from 


{llur 1 M - “7 [i]} {a} = o 




(13) 


Thus, for the k, m-th natural frequency, the eigenvectors give the displace- 
ment components as 

N -i 

U k = B 2 k X 2 + Bi k (X - 1) + B mOk + B mnk cos j c ° sm ^ 


N 


V k = I C mnk sin sinm0 

n = 1 ' ' 


w k = Z Amnk sin cos m0 

n = 1 v * ' 
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Generalized Added Mass of Liquid 

For radial shell vibrations that are proportional to cos m0» the liquid 
exerts a radial pressure which may be regarded as an apparent mass M mn t n 
which is added to the cylinder mass. We now derive an expression for this 
quantity. 

The fluid is assumed to be nonviscous, ir rotational, but compressible. 
The velocity potential corresponding to the mn-th component of shell motion 
for small disturbances is governed by the wave equation which for periodic 
motion is 


V^ < j > mn + ^hnn - 0 
c 0 


The boundary condition on the wall is 
9 < l ) mn _ 9 w mn 


at r = a (w mn outward positive) 


(14) 


9r 9t 

where the mn-th component of radial shell motion is 

w mn = A mn (t)f n (x)cos (m0) (15) 

and we assume that f n (x) is an orthogonal set of functions in the interval 
[O, i] . The boundary condition at the bottom is; 


ifeE=0 at x = 0 (16) 

9x 

When frequency of excitation is much higher than the leading few liquid 
surface sloshing frequencies in consideration, the free surface condition 
can be approximated by 

= 0 at x = h (17) 

By separation of variables, a particular solution can be constructed to 
satisfy the boundary conditions on the wetted surface [Eqs. (14)and(l6)] as; 
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*1 = Dkn^mk cos (^)cos(m0) ( a — - 1 ) 

with Fourier expansion of f n (x) being 

*„<*>= I Dkn«* (^r) 

k * 0 \ n / 

where 



and 


^mk 


Im ( Ck a) 


k > 


loh 

itCq 


Jm(hf) 

mk SkJm(lk) 


for 


k < 


coh 



Next consider a complementary solution, satisfying homogeneous boundary 
conditions on the wetted surfaces so that the net velocity potential satisfies 
the approximate free-surface condition: 
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where 


and 


is the j-th root of Jjn^mj) = 0 
C j (“> w ) = cosh (n m j for jirnj >“ 

c j (a’ W ) = COS (^J a) f ° r < 


where 


f"l 2 u>2 a 2 

r|m ) = LI 1 * 1 " 3 ’ ~cT. 


E kj = 




%0=2/ f R Ok(j.w)d(j) 


oo 


B 


X (“1)^ f^kn^kj 
k = 0 


mjn 




OO 


£ ( “1)^ Dkn^kO 
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The total velocity potential is given by 
<j>mn ~ <f>l + 4*2 

and the linearized pressure loading, given by the Bernoulli equation, is 

^4mn 


Pmn “P 


3t 


r = a 


As a result, the n' -th component of generalized force p mn ' n is given by 
integration of the loading with a weighting function f n i(x): 


i 


/ f n'( x ) Pmn^ 


Pmn' n 


X. 

/ 


’^mn'nP g h s c °s md 


d 2 A 


mn 


f 2 i (x) dx 


dt 


where the generalized apparent mass is 

oo 

no 1 

M mn *Tl “ 


Rmk(l) w ) DknDkn 1 


r 

I 

Ps h s Lk = 1 

, ? E kODkn <-U k £n'50m/ cos (^)] + (jT c0g «h 

K “ * ^ r* 


J 0^0> - 

Don' 


loJido) 


[' 


D0n 6 m0 + ( 1 " 5 m0) Do n Don’ * 


oo 


- X J m(Pmj) c nj B mjn 

j=l 


} 


in which 



r] 

(18a) 
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f n i(x) dx 



/ cos (|f) f n'< x ) dx 

« 0 A* 

E n , - - = DQn' 

/ fn'(x) dx 

0 

Note then that the mn-th component of pressure loading on the shell is 


2? V d 2 A mn 

Pmn ~ ^ Pmn'n^n'( x ) ? ^ p s h s M mn i n cos (m0 )f n t(x) — 

n' = 1 n' = 1 dt^ 

(18b) 

Thus, M mn i n is the coefficient of p s h s (d 2 /dt 2 )(A mn ) in the n'-th component 
of shell motion exhibiting a cos (m0) mode of vibration. To use this generalized 
coefficient in Equation (11), we note from Equation (2c) that f n '(x) assumes the 
form 



Forced Axisymmetric Response 

Now, assume the tank is excited axially. A similar Fourier process 
will be utilized for the solution of this case, so that the same shell displace- 
ment forms [Eqs. (2)] are again used, except that a superscript p will be 
incorporated to designate the forced motion. Further, the response will 
occur at the excitation frequency ft, rather than at a natural frequency ft^. 
Finally, in this section we dispuss only axisymmetric motion (m- 0), which 
is of the most practical concern in a linear formulation. 
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The boundary conditions on axial displacement become 
At X = 0 , U = -X Q cos (tot) 

At X = - , 3It = M**£2 2 U 

a dX 

To satisfy Equation (19), one has, from Equation (2a): 



(19) 

( 20 ) 


Additional equations corresponding to Equations (7a), (8), (9), and 
(10) can now be written as before. For brevity we do not write these 
equations since they are identical to those mentioned except that the 
coefficients carry the superscript p and £2 replaces £2^. Finally, an 
equation corresponding to (11) can also be written for the forced response. 
However, we must first discuss the difference in pressure loading that 
will occur. 


As a result of the axial motion of the tank bottom, the liquid exerts 

B 

a generalized radial pressure loading pQ n on the tank wall in addition to 
the apparent mass pressure p^ in . That is, for forced motion the total 
pressure can be expressed as 


p p = p p + p* (21) 

*r man ^0n 

The apparent mass pressure pP^ has already been determined by Equation (18b) 
providing that A^ n is utilized in lieu of A mn . The additional pressure loading 
remains to be determined. This pressure can be derived from a one- 
dimensional velocity potential which satisfies the wave equation along 


with the boundary conditions: 
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4> B - 0 at x - h 


— * xm sin cot at x = 0 
9x 0 

The velocity potential which satisfies these conditions can be found as 


insinutsin [ — - (h - x)l 

l c Q J_ 

go / wh \ 

— cos { — ] 

c o V c o/ 


Thus, the generalized loading (which results from the application of the 
Fourier process) relative to the n’-th tank-displacement component can be 
expressed as : 


P R,.- JB- 


g m pw^xq cos (wt) sin (h - x)^ f n i(X) 

On _ 1 / V c 0 L 


= 4 / 


oo / coh \ 

— cos — 

c 0 \ c 0 / 


where 


The amplitude of this loading becomes 




(™h) 


where we define 


= -JL J s i n T— (H - X)1 sin(\ n , 

a 2 J 0 Lc 0 J 


X) dX 
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We now are in position to formulate the last equation required for 
the solution of the forced response. Combining Equations (18b), (21), and 
(22), along with Equations (2) with superscripts p on the coefficients, into 
Equation (1c), and collecting the coefficients of the sin (n ! 7rx/i ) terms, 
there results 


oo oo 

v + B^ 0n , - £ VO n .J +rn £ C^ n 6„, n 

L n = 1 -* n = 1 


oo 


H 


s 2 , 2 >2 . 1 - v‘ 


+ I A mn8n , ni 1 + “72 (Kn + m } + 

n = 1 


H c 


"x 2 (E± . ¥g ^ 

L n \ 2E 2Tra 2 E j 


+ m 


2 p 0 


EJ H 


2 oo 

+ 1 m2 Z £ ^ x 3n'n A n 


n = 1 


mn 


(i - v 2 )n 2 I 

n = 1 


oo 


^n'n + ¥mn'n 


= (l-v 2 )pf nI (23) 


The equations for determining the forced response can now be put in the 
following matrix form 



(24) 


2 = (1 - v 2 ) £2 2 


where 
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and the elements of the 


matrix are 


F rn' = (1 - v 2 ) P^ n , 

The lull and IrI matrices are the same as in Equation (12). Here, 
however, the column matrices of response coefficients are different 
from those of the natural modes. 

For given values of Xq and , one can solve for the forced response 
by the following inversion: 


(a p } 

L mnJ 

(B p l 
\ mnj 

{ c iU 


rP 

c m0 


• = X 0 flul - 5 2 iRlj _1 


{ F r} 

{»} 

{“} 


Q 
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THEORETICAL AND EXPERIMENTAL RESULTS 

Figure 2 shows the variation of natural frequency with liquid depth for 
several modes of a steel tank whose geometry is given in Figure 1. Each of 
the three parts of Figure 2 is for a different top mass. Theoretical computa- 
tions were performed for the largest mass (Fig. 2a) only. All symmetric 
modes in the given frequency range were obtained, although only one non- 
symmetric mode is indicated. Of course, many other nonsymmetric natural 
modes existed in this range (note that experimental data for two nonsymmetric 
modes are given in Figures 2b, c). 

Natural frequencies for the symmetric modes were determined ex- 
perimentally by detecting peaks in the pressure at the center of the tank 
bottom and/or detecting peaks in the output acceleration of the top mass. 
Although Figure 2 indicates a condition of zero ullage pressure, up to pg = 

10 psig had to be used in order to obtain data for the symmetric modes. This 
was necessary to prevent the occurrence of instability in some nonsymmetric 
mode, in which a simple linear symmetric response no longer was present. 
Fortunately, this procedure was possible since the frequencies of symmetric 
modes were determined to be independent of ullage pressure. We emphasize, 
however, that frequencies of nonsymmetric modes are highly dependent on 
ullage pressure, so that data for the k = 1, m = 10 mode (as well as the k = 1 , 
m = 1 3 mode) in Figure 2 were taken at p 0 = 0. Data for these modes were 











b. W 0 = 22.32 lb. 

Figure 2. Natural Frequencies Of Partially Filled Tank With Top Mass 
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Natural Frequency, cps 
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taken as peaks in the wall response at the antinode of the axial wave form. 

Of course, the position of this antinode shifted with different liquid depths. 

Theoretical data for Figure 2a were obtained by the use of a five- 
term (n = 1 to 5) expansion in Equation (13) and computing the eigenvalue 
for the resulting 18 X 18 matrix. Some deviation between theory and exper- 
iment can be seen to exist in Figure 2a. However, the use of more terms 
in the expansions would reduce this deviation. This can be seen from the 

three points computed for n = 10 terms for the first three symmetric modes 
at a depth of h/f = 0. 55. Since the five -term expansions appeared to give a 
sufficiently good comparions between theoretical and experimental results, 
most of the computations were thereby limited in order to reduce the 
required digital computer time. 

Several interesting observations can be made from the data in 
Figure 2. Although the frequency of the m = 10 mode is considerably below 
those of the symmetric modes throughout most of the depth range, pres- 
surizing the tank can raise the nonsymmetric mode above the lowest sym- 
metric mode. For low liquid depths, it can be seen that the first symmetric 
mode represents the first coupled axial top mass -shell mode with only small 
liquid effects, while the second symmetric mode represents the first coupled 
liquid-shell mode with only small top-mass effects. However, these roles 
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of the first two modes are interchanged for greater liquid depths as indicated 
by the dashed lines in Figure 2a which represent the respective decoupled 
modes. 

From Figure 2 it can be seen that variation in top mass had no effect 
on the nonsymmetric mode but a strong effect on the first two symmetric 
modes for a middepth range of liquid. This is shown further in Figure 3. 
Finally, the influence of ullage pressure on frequencies of the nonsymmetric 
linear mode was not measured since this has been determined by various 
previous investigations. 

Figure 4 shows a comparison of theoretical and experimental axi- 
symmetric forced response for a frequency range which includes the first 
two modes with a liquid depth of h/f = 0, 69. Although the data were taken at 
Pq = 1 0 psig, the results are independent of pressure, as has already been 
mentioned. Here, the acceleration amplification of the top mass was chosen 
as a comparison parameter, although the liquid pressure at some point in the 
tank could have been used just as well. An intermediate liquid depth was 
chosen as a worst possible condition for using a given number of series terms 
in the theoretical computations. That is, previous work indicates that the 
most serious distortions of tank axial mode shape from a half-sine wave 
occurs at intermediate depth ranges. 

Theoretical points were determined from Equation (25) by 
means of a digital computer. Of course, the net sum of the forced 




Figure 3. Effect Of Top Mass On Natural Frequencies 
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Figure 4. Response Of Top Mass Relative To Input Acceleration For Axisymmetric Wall Response 
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axial displacement components u at x = 2 form the part of the theoretical 
solution which is used in Figure 4, Basically, five terms were 
used, although ten terms were used for part of the range as indicated. It 
can be seen that more terms reduce the discrepancy between theory and 
experiment. Also, the location of the theoretical and experimental resonance 
points corresponds with the location of the respective natural frequencies 
for the first two symmetric modes in Figure 2a. Thus, the agreement 
between theory and experiment can be made as good as is desired, and the 
most terms in the expansions are required at intermediate depths. 

A final comment may be made with regard to the use of polynomial 
functions in conjunction with a Fourier series (see Eq. 2a). It is recognized 
that convergence of the solution may be of some concern. However, the 
relatively good comparison between experimental and theoretical results 
tends to refute the latter concern, at least in the range of parameters 
investigated. 
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APPENDIX 


Elements of Submatrices in Equation (12) 


n' designates row and n designates column of the respective elements 


(\£ + m 2 ) 2 6 n « n 


U ln r n ~ ^ + 19 (^n + m ) 


? 2 

1 - v 2 po/^n , ?\ , ^g^n , , pga ?v 

U7 _ [eVT m ) n ' n ' TTzi 6 " n + "I" m X 3n'nj 


^ln'n “ ”^n v ^n'n » ^ln'n “ m ^n'n 


Yin' = v X 0 n' » ^i n t - vXln' , R ln'n - ^n'n + Mmn'n 


^2n'n “ " v ^n^n'n » ^2n’n - ^n + ^ 2 ^n'n 


^2n'n " “ m ^-n^n'n 


(1 + v) 


, Y 2n t = m 2 Xi 


**!»• ^ 


,1 = m 2 X ?ri i - 1 ■ ; V - , S, n . 


2 n r n °irn 


F > 2n' = *ln' » ^2n' “ 


U3 n 'n - m ^n'n » v 3n'n " _ m ^n^n'n 


(1 + v) 


, n = m 


2^2 


(1 - v)' 


n 2 


Aa'n » ^3n' “ m * 0 n' 


(1 + v) 


Z 3n ' = mX] n i — » T 3 n 'n = Vn 



